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Abstract: We introduce generalized energies for a class of Uq{Dn^) crystals 
by using the piecewise linear functions that are building blocks of the combina- 
torial -R. They include the conventional energy in the theory of afSne crystals 
as a special case. It is shown that the generalized energies count the particles 
and anti-particles in a quadrant of the two dimensional lattice generated by 
time evolutions of an integrable Di^' cellular automaton. Explicit formulas 
are conjectured for some of them in the form of ultradiscrete tau functions. 

1. Introduction 

Let Bi be the crystal of the Z-fold symmetric tensor representation of the quantum 
affine algebra Uq{Dn^) [HJIU]. The combinatorial R : a; (g) ?/ M- y' (g) a:' is the isomorphism 
of crystals Bi g) B^ — > B^ (g) Bi corresponding to the quantum R a.t q = [T^. In 
Ref. [15) . an explicit formula of the combinatorial R was obtained in terms of several 
piecewise linear functions gi{x ® y) G Z on Bi (g Bm- See Theorem 12.11 Among them 
is the local energy, which plays an essential role in the theory of affine crystals |10) . 
The family of piecewise linear functions {^i}, which we call generalized local energies 
in this paper, are ultradiscretization of the subtraction-free rational functions that have 
emerged as building blocks of the tropical R [TSl [TB] of the geometric crystal |3] . They 
may be viewed as local energies in a principal picture rather than in the conventional 
homogeneous picture. 

From the local energy, one can form the integer-valued function called energy on the 
tensor product V = Bi^ g) • • • g) B/^ . Its generating function is the one dimensional 
configuration sum that originates in the corner transfer matrix method [1] [2] . 

In this paper we introduce generalized energies £g. : V ^ Z>o corresponding to g^'s, 
and study them from the viewpoint of the integrable cellular automaton of type Dn"^ 
[51 [7]. The latter is an integrable Uq{Dn'^) vertex model at q = 0. It is a dynamical 
system on V equipped with commuting time evolutions {T/}/>i. Elements of V are 
naturally regarded as arrays of particles and ant i- particles, and Ti induces their factorized 
scattering involving pair creation and annihilation. See Examples 13.11 and 13.21 

Our main result is Theorem l4.11 which states that . (p) = pg^ {p) for any p = pi g) • • • g) 
PL ■ Here pg. [p) is a counting function giving the number of certain particles and anti- 
particles specified by gi in the region p.3|) under the time evolutions p, Too{p), T^{p), . . .. 
As such, the counting functions are non-local variables attached to a quadrant of the 2 
dimensional lattice. However, it will also be shown in Theorem 13.51 that the combined 
data {pgiipi g) ■ • • g) Pj) \ j = k — l,k} in turn reproduces the local variable pk S -B/^. 
completely in agreement with the spirit of the corner transfer matrix method. Therefore 
the joint spectrum {Sg^ipi g) • • • g> Pfc)} of the generalized energies with 1 < < L is 
equivalent to p — pi (E) ■■■ (S pl G 'P itself. This extends a similar result on type An ^ 
(Proposition 4.6 in Ref. [T^) which is related to the katabolism (35]. A supplementary 
result (Proposition 1331) is parallel with Theorem 14.11 and treats generalized local energies 
with opposite chirality (cf. Remark l2.5p . 

The layout of the paper is as follows. In Section [2l generalized (local) energies are 
extracted from the piecewise linear formula of the combinatorial R [15] . In Section O 
the integrable Dn^ cellular automaton |6l [7] is recalled and the counting functions are 
defined. In SectionUl the main Theorem 14. II of the paper is stated and proved. In Section 
[5l aspects related to combinatorial Bethe ansatz are discussed. In Section [5711 we give 
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the inverse scattering formahsm of the Dii ' cellular automaton like Ref. jMj . In Section 
15. 2[ we conjecture piecewise linear formulas for some generalized energies in terms of 
ultradiscrete tau functions. This is also motivated by the ^4^""^'' case [17j . where analogous 
results have led to a piecewise linear formula for the Kerov-Kirillov-Reshetikhin map 
[IB] . Although the conjecture is yet to cover the full family of generalized energies, 
the last one (|5.8p is already rather intriguing. We expect that the extention and the 
solution of Conjecture 15.31 will uncover an interplay among combinatorial Bethe ansatz, 
ultradiscretization of the DKP hierarchy 8: and the bilinearlization of the tropical R 

m- 

2. Generalized energies for dI^ '' crystal 

2.1. Crystals and combinatorial R. Let us recall the basic facts on crystal and com- 
binatorial R briefly. For a more information, see Refs. [Tni[TTl|2] and [12]. For a positive 
integer I, let 

n 

(2.1) Bi = {C = (Ci, . . . ,Cn,C„, • • • ,Ci) e zf„ I + Q = I, CnC„ = 0} 

i=l 

be the crystal of the Mold symmetric tensor representation of Uq{Dn ^) [9 . We assume 
n > 3. As for the functions e^, ipi, the tensor product rule and the action of Kashiwara 
operators Ci and /i (0 < i < n) , see Ref. [19] . 

The affinization of the crystal Bi is defined by AfF(i?;) — {b[d] | d G Z, & G Bi} with 
the crystal structure ei{b[d]) = {eib)[d + 6io] and fi{b[d]) — {fib)[d — Sio]- We call b and 
d the classical and the afSne part of b[d], respectively. There exists the unique bijection 
(crystal isomorphism) Bi®Bm -> B^^Bi that commutes with all Kashiwara operators. 
It is lifted up to a map AE{Bi) Aff(i?„i) ^ Aff(_Bm) Cx) AfF(_B;) called the combinatorial 
R, which has the following form: 

R : AS{Bi) ® AS{Bm) Aff (B„0 ® Aff (5,) 

b[d](E)b'[d'] ^ b'[d' + H{b(E)b')](E>b[d-H{b(E)b')], 

where 6 6' ^ b' ^ b under the isomorphism Bi B^ ^ Bm (K> -B;. Q The quantity 
H[b (E) b') is called the local energy and determined up to a global additive constant by 

r H{b ® fe') + 1 if « = 0, Mb) > £oib'), Mb') > eo{b), 
Hii,ib®b'))= I H{b(^b')-1 ifi^O, Mb) <So{b'),Mb') <Soib), 
[ H{b®b') otherwise. 

The Yang-Baxter equation 

(2.2) (i? ® 1)(1 ® R)(R ®\) = {l® R){R ® 1)(1 ® R) 
is satisfied on Aff(B() Aff(B„) ®AS{Bk). 

2.2. Generalized local energies. Let us give an explicit piecewise linear formula of the 
combinatorial R that originates in the tropical R for geometric crystals of type Dn^ |15) . 
First we make a slight variable change. The set Bi (|2.1I) is in one to one correspondence 
with another set 

B'l — {x = {Xi^ . . . ,Xn,Xn-l, ■ ■ ■ ,Xl) G 7/^"^^ \ Xi,Xi > for 1 < i < 71 — 1 , 

71-1 

(2.3) Xn > -min(a;„_i,x„_i), ^(xi + Xi) + x„ = 1} 

i=l 



This classical part of the combinatorial R will also be referred as combinatorial R and denoted by 
R{b 6') = b' ig) b. 
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by the relations 

(2.4) X, = C^, x^=Z, (l<j<n-2), 

(2-5) Xn — l = Cn— 1 Cni •^n Cn ^ Cm -^ri — 1 Cri — 1 ^" Cm 

(2.6) Cn = max(0,Xn), C„ = max(0, -a;„), 

(2.7) Cn-i = a;„_i +mm(0,x„), C„-i =Sn-i +min(0,a;„). 

Note that be negative. We naturally use the notations like x[d\ £ AS {B'l) and 

R{x[d] ® x'[d']) = x'[d' + H{x a;')] ® x[d - H{x (g) a;')], etc. Set 



(2.8) ^(C)=^(C. + C.) (CeBO, e{x) = J2ix^+x,) + Xr. (X6B0, 

1=1 i=l 

SO that ^(C) = l{x) = 1 ioT C € Bi and x e B'l. 

Let X = (xi, . . . G B'l and y = (yi, . . . ,yi) G B^. On the pair (a;, y) we introduce 
mutually commuting involutions cti , cr„ and * by 

(2.9) ix,yr = {x'^\yn. {x,yr- = {x''-,y^-), {x^y)* = [y* ,x*), 



CTl 

(2.10) * 



i > Xl, 

a;i ■< — > Xi (1 < j < ft — 1). 



The coordinates not included in the above rules are left unchanged. These involutions 
are naturally defined on (^, C) & Bi x Bm as well by the correspondence (|2.4l) - (|2.7p . 
For instance, oiie has (|,C)* = {CO with C* = (Ci, • ■ • , Cn-i, Cn, C„, Cn-i, • ■ • , Ci) for 

C = (Cl, • ■ • : Cn, Cn: • ■ • , Cl)- 

For any function g = g{x,y), we write g"^^ = g"^(x,y) = g{x"^ ,y"^), etc. Introduce 
the piecewise linear functions Vi — Vi{x,y) and Wi = Wi{x,y) for < i < n — 1 as 
follows. 

(2.11) V, = max {{0^,J,eyl < j < n - 2} U H.j,vlj\'^ < J < n}) , 

(2.12) Wo = 2Vo, Wi^Vo + Vo^\ = K-i + K-i, 

(2.13) Wi = max {V^ + V*_^ - y,, V,-i + V* - x,) + min(x,,|/J, 

where 2<z<n — 2in the last line. The functions Oi^j — 9ij {x, y), 9[ ^ ~ 6[ ^ (a;, y), rjij = 
Vij{x,y),vlj = Vt,j{x,y) are defined by 



^{x) + ^ {Vk- Xk) for 1 < j < i, 
k=j+i 
j 

Ky) + J2 (^fc -Vk) ioT i + l <j <n-2, 

^ k=i+l 

i 3 

%^j{x, y) = l{x) + ^(yfe - Xk) + ^{yk -Xk) for 1 < j < n - 2, 



fc=i fe=i 



Table 1 . Transformation by cri , cr„ , * and R. 





Vo 


V{1 <i<n~2) 


Vn-l 


W^{l<i<n-l) 


* 

R 


Vo 
Va 
Vo 


V 
V 
V* 
W, - V* 


14-1 





l{x) + {Vk- ^k) + Vj - Xj for 1 < j < i, 

k=j+i 





Vi,3 (x, y) = \ (-{y) + ^ {xk-yk) + yj-xj for i + 1 < j < n - 1, 

k=i+l 
n-1 

^(y) + ^ {xk- Vk) + Xn for j = n, 

^ fc=i+l 

^{x) + ^(ffc - Xk) + ^{yk - Xk) + Xj - Ijj for 1 < J < n - 1, 



vlj {x, y) 



k=l 



k=l 



^{x) + 5„,-i {e{x) - e{y)) + J^ivk - xk) + J2iyk- xk) 

k=l 

for j = n. 



k=l 



Theorem 2.1 (Ref. [TS], Theorem 4.28 and Remark 4.29). The image y'(g)x' = R{x(E)y) 
of the combinatorial R is given by 

x',^x, + V^^~V*, x[^x, + VU + W,-V* ^W,-i (l<i<n-l), 

(2.14) < = 2;„ + v;*-i -14-1, = yn + K-i -v;*_i, 

y[^y^ + V,^l+W,-V,-W,^l, y[=y^ + V,.i-V, {l<l<n-l). 

Moreover, the local energy is given by 

(2.15) H{x®y) = Vo{x,y) 
up to a constant shift. 

The functions Vi , . . . , Ki-i , Wi , . . . , W4-i and ui , cr„ and * of them are relatives of the 
local energy. In addition to the involutions cri,CT„ and *, the combinatorial R naturally 
acts on them by {RVo){x, y) — Vo{y' , x') with y'(^x' — R{x®y), etc. Their transformation 
properties under ai,an,* and R are summarized in Table IT] |15j. These involutions are 
commutative, thus for instance R{Vq^) = (i?(Vb))'^^ = Vq^. 
Due to these properties, there are a few simplifications in (I2.14p as 

x[ = xi + Vo-V{, x[^xi + Va^' -V{, 

(2.16) y[ = yi + Vo^'-Vi, y[^y, + Vo~Vi. 
We write 

(2.17) ui^{l,0,...,0)eBi. 
By using Theorem 12. 11 one can show for any ^ G Bm that 

(2.18) Bi (X) Brn 3 Ul(g)C ^^ Urn ® S.' ^ Bm ® Bi li I > m 

for some ^' under the combinatorial R. In particular 

Ul^Um—Um® Ul 



5 



holds. The functions in Table [T] attain their maximum Vq ~ V^^ = Vi = V* = l + m and 
Wi = 2{l + m) tor 1 < i < n — 1 at {x, y) = {ui,Um)- 

For (E) ( ^ Bi (g) Bm , let a; G B[ and y G be the elements corresponding to ^ and 
C, respectively. We set 

(2.19) ^'.(C®C) = m+i{0-V,{x,y) (0 < i < n - 1) 

(2.20) vl-{S,®0 = m+m-V^'{x,y), 

(2.21) <(c®c) = m+m~v:{x,y) (i<z<n-i), 

(2.22) w,(C®C) = 2i{^) + 2t{0-W^{x,y) {l<i<n-l), 

and call them generalized local energies. Note that — Vn-i — They are 

building blocks of the piecewise linear formula of the combinatorial R (|2.14p . From the 
above remark, generalized local energies are all nonnegative and normalized so that 

(2.23) g{ui (g) Um) = for any g — Vi, , v* and Wi. 
For C = (Ci, • • • , Cn, Cn, ■■■Xi) ^ Bi, we introduce 



(2.24) 


a(C) 


= C2 + 


•••C« + C„ + --- + C2 + 2Ci - 


:^(C)+Cl-Cl, 


(2.25) 




= C2 + 


•••C«+C„ + --- + Ca+1+Cl 


(0 < a < n - 2), 


(2.26) 




= C2 + 


•■•+Cn-1+Cn+Cl, 




(2.27) 


7<_,(C) 


= C2 + 


•■•+Cn-1+C„ + Cl, 




(2.28) 


7«>a-«a(C) 


= C2 + 


•■• + Ca + Ci (l<a<n-2) 


1 


(2.29) 


7^„'i (C) 


= 0. 







Note that jy„{C) = a(C). 

Lemma 2.2. Let ^ = (Ci, ■••,?!) e Bi and C = (Ci, •■•,(!) e B,„. Set (' ® ^' = 
R{S, (S) C) ^ <8) Bi. For (hence I as well) sufficiently large, the relation 

g(C®C) -7^(0 + a(C') -7^(0 

is valid for g appearing in h2.25\) - P2.29\) . where the left hand side with g — Wa — Va is to 
be understood as Wa{^ ® C) ~ ''^aCC ® 0- 

Proof. One can check that > m — Ci + is sufficient to guarantee that Vq ()2.1ip is 
equal to r?Q ^ = ^ + Ci ~ Ci- This implies that wo(C <8> C) — rn — (i + Ci = 7i'o(C) showing 
the g = Vq case. All the other cases are deduced from this, (|2.14p and (|2.4p - (|2.7p without 
using the concrete forms of 9ij, 9[ j, rjij and r]^ j. □ 

For g = Wa {1 < a < n — 2) , a.n analogue of Lemma holds with 

wa{^<»0 = 7»„(C)+2a(C')-7»JC'), 

IwAO = 7t«a-fa(C) +7i;a(C) 

(2.30) = a(C)+C2 + --- + Ca-(Ca + --- + C2)- 

As ^1 gets large, (' stabilizes since it is a piecewise linear function of staying in 
a finite set B^. {£,1 > m seems sufficient for the convergence.) Therefore Lemma [2.21 
ensures that all the generalized local energies g{^®C) are well defined in the limit — > oo. 

2.3. Generalized energies. For p = pi ® ■ • ■ ® pl G Bi-^ ® ■ ■ ■ ® Bi^ , define pj'-* G Bi. 
{i < j) by 



(Bi, ■ • • ® Bi^_^ ) ® Bi^ ^ Bi^ ® {Bi^ ® • • • ® Bi^_^ ) 

f ®p',®---®p'^_i, 



(2.31) Pi® ■ ■ ■ ®Pj-i(S)Pj ^ pf ®p',®---®p' 
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sending pj to the left by successive applications of the combinatorial R. We set pj" — pj. 
For any generalized local energy in (|2.19|) - p.22p . we define the generalized energy of 
p ^ Pi ® ■ ■ ■ ® PL & Bi^ ® ■ ■ ■ ® Bi^ by 

(2.32) £,{p)= fffe^^r'^) 

0<i<j<L 

by taking po — ui with sufficiently large I. This is well defined (finite) due to Lemma [2?2] 
and the comment following it. In the rest of the paper, we will simply write pq — Uqo G 

Boo- 

When g — vo, p.32p is the energy introduced in Refs. [18] and up to a sign and a 
constant shift. If furthermore h, . . . ,Il are all equal, then pj*^^'' = Pi+i holds and ()2.32p 
reduces to 

^vo (P) = (-^ ~ *) ^0 {Pi «) Pi+l) ■ 

0<i<L 

Its generating function q^^o (p) is a version of the one dimensional configuration sum 
going back to Refs. ^ and pj, which is the essential ingredient in the corner transfer 
matrix method. 

Any quantity G{pa(S)Pa+i(S)- ■ •®pp) will be said R-invariant if G(- • ■®pi®pi+i®- • • ) = 
G(- • • ® R{pi ® Pi+i) ® • • • ) for any a <i < {5. 

Remark 2.3. Due to the transformation property under R in Table [1] the generalized 
energy £g{p) is i?-invariant for g = wq, Wn-i, '^'^^-i and wi,.. . , On the other 

hand, £g with g = ui, . . . , u„_2 and vl, . . . , w*_2 are not i?-invariant. 

Let us depict the relation R{b ® c) = c®b a,s 



b c 




c b 



Then the Yang-Baxter equation p.2p takes the well known form: 




The defining relation (|2.3ip of looks as 

Pj-l Pj 

1 

Pi 

(2.33) 

Remember that each vertex is associated with various generalized local energies g(b®c). 
Let 

(2.34) ig = lg{p^ (g) ■ ■ . ^ Pj ) = Y 9{Pk ® Pf^^^) 

i<k<j 

be the sum of generalized local energy g over all the vertices in p.33p . Then the gener- 
alized energy (I2.32p is expressed as 

Sgipi ® ■ ■ ■ ® Pl) = £g{pi ® ■ ■ ■ ® PL-i) + IgiUao ®Pl ® ■ ■ ■ ® Pl) 

(2.35) = Y, Ig{Uoo ®Pl ® ■ ■ ■ ®Pj-l ® Pj)- 

1<3<L 
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From (I^JS)) and it follows that 

(2.36) Ig{Uoo Uoo (E)P1 ■ ■ ■ 'S) Pj) = Ig{Uoo Pi (E> ■ ■ ■ ®Pj), 

(2.37) £g{Uoo<^Pl<S)---®PL) = £g{Pl<S)---®PL)- 

Lemma 2.4. /n \2.33\) . the following quantities are R-invariant as the functions of 
Pi® • ■ • ®Pj-i- (i) The element p^j'K (ii) Ig {2. 34]) for any g = Va {0 < a < n ~ 1) , Wa — 
Va ^ a, < n — 2) , and Vq^ . 

Proof, (i) This is due to the classical part of the Yang-Baxter equation ()2.2p . (ii) The 
i?-invariance of I^g follows from ()2.19p . ()2.15p and the affine part of the Yang-Baxter 
equation. Let y,y' G B[, be the elements corresponding to Pj,p^j ^ G Bi^, respectively. 
From (|2.16p . we have y[ —yi= Ivi — ho- Since the left hand side is i?-invariant by (i), 
this relation implies the i?-invariance of . By similarly using the i?-invariance of % — 
and y'^ — y-i in (j2.16p and (|2.14l) , one can verify the i?-invariance of the other Ig . □ 

We note that Lemma [2.41 is applicable to the situation i — i.e., pi ® ■ ■ ■ ® Pj-i ~ 

Uoo ® Pi ® ■ ■ ■ ®Pj-l- 

Remark 2.5. Lemma HH] (ii) does not concern 1^^, . . . , _j • In fact the proof does not 
persist since V{, . . . , V*_2 are not contained in any difference of the components of y' 
and y in ()2.14p . Similarly, Vi, . . . , Vn~2 do not appear in the differences of x and x' . This 
chirality of the combinatorial i? is a characteristic feature of the Dn'^ case. In contrast, 
Ig's with all the generalized local energies g for An"^ (so called ith (un) winding number 
[T7] ) are i?- invariant. We shall come back to this point again in Section 

The following proposition and its proof are parallel with Lemma 4.4 in Ref. for 
type aI^^ . 

Proposition 2.6. For g = Vi {0 < i < n — l),Wi {1 < i < n — 1), v^^ and v^_i, the 
generalized energy £g{pi ® ■ • ■ ® pl) 32\) is equal to the sum of the generalized local 
energy g attached to all the vertices in the following diagram (L — 3 example): 



Pl P2 P3 




Proof. We invoke the induction on L. For L = I, one has £g{pi) = g{uoo ®Pi), and 
the assertion is obviously true. We illustrate the induction step from L = 2 to L = 3. 
Consider the following identity obtained by successive applications of the Yang-Baxter 
equation: 

Uoo Pl P2 P3 Uoc Pl P2 P3 





Here •.Ci^di stand for the values of g at the attached vertices. By the induction as- 
sumption, the sum of the three • in the left hand side is equal to £g{pi ®P2)- In view 
of the recursion relation (12.351) . we are to verify ei + 62 + £3 = Ig{uoo ® Pi ® P2 ®P3)- 
By the definition, Ig{uoo ® Pi ® P2 ® Ps) ~ di + d2 + d^ in the right diagram, where 

di = g{p2 ®p'3^)i d2 = g{pi ®p'3^),d3 — g{uoo (81 P3^'). Thanks to the i?-invariance of Ig 
in Lemma 12.41 (ii), this is equal to ei -I- 62 -f 63. □ 
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3. Integrable d';P cellular automaton 

3.1. States and time evolution. Let us recall the integrable D^P cellular automaton 
associated with i?; O [7]. Consider the crystal Bi^ • • • ® Bi^. Its elements are called 
states. We regard each component {(^i, . . . Xi) G Bi as a capacity / box containing Ca 
particles a and anti-particles a for 2 < a < n, and furthermore Ci extra I's which we 
call bound pairs. The remaining I's represent empty space in the box. Thus ui stands 
for an empty box. The indices a, a will be referred as color of particles and anti-particles, 
respectively. A state pi^ - ■ ■ ®pl G Bi^ ® ■ ■ -^Bi^ represents a configuration of particles, 
anti-particles and bound pairs in an array of boxes with capacity li, . . . ,1^. Because of 
the constraint CnC™ =0 p.ip . particles n and anti-particles n do not coexist within 
a box. We shall denote the element (3,0,1,0,2,0,1,0) £ Bj of d'^\ for example, by 
1113442, etc. 

For a positive integer I, we define the time evolution Ti{p) = p'^ ® ■ ■ ■ ® p'j^ of a state 
p = Pi® ■ ■ ■ ®PL by 

Ui ® pi ® ■ ■ ■ ® PL — P'l ® ■ ■ ■ ® p'l ® i 

under the isomorphism Bi (g) {Bi^ ® ■ ■ ■ ® Bi^) ~ (i?;^ • • • Cg) Bi^) (g) Bi. Here ^ G -B/ as 
well as Ti (p) are uniquely determined from p by the combinatorial R. It can be shown 
that 

(3.1) £, = ui a pj = ui- for L' < j < L with sufficiently large L — L' . 

The time evolutions {T;} form a commuting family, and T; stabilizes as I gets large, 
which will be denoted by Too. 

When li — ■ ■ • = II ^ 1, T^c is factorized as 

with Ka given by the following algorithm (we understand a = a) [7] . 

(1) Replace each 1 by a pair a, a within a box. 

(2) Pick the leftmost a (if any) and move it to the nearest right box which is empty 
or containing just a. (Boxes involving the pair a, a are prohibited as the desti- 
nation.) 

(3) Repeat (2) for those a's that are not yet moved until all of a's are moved once. 

(4) Replace the pair a, a within a box (if any) by 1. 

In the above, taking some 6 (^ 1) g Bi away from a box means the change of the local 
state 6 — 1. Similarly, putting 6 (7^ 1) e Bi into an empty box means the change 1 — > 6. 
The steps (1) and (4) can be viewed as pair creation and annihilation, respectively. 

Example 3.1. We consider d'^'' and states from _B®^ with L = 48. Successive time 
evolutions of the initial state on the first line under Too is presented downward, where 
61 (g) ■ • • (g) 6l is simply denoted by a horizontal array 6162 ... fo^. It shows a collision of 
solitons with amplitudes 6 and 3. 

111334322111112241111111111111111111111111111111 
1111111 1 13343221 12241 1 11111111111111111111111111 
11111111111111133 432 T2 4 1111111111111111111111111 
111111111111111111111 34TT4 44 1 1111111111111111111 
1111111111111111111111111422223 4 4411111111111111 
1111111111111111111111111111422111223 44411111111 
1111111111111111111111111111111422111111 223444 1 1 

For general. Too still admits a similar, although slightly more involved, al- 

gorithm. We omit it here and give an example instead. 
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Example 3.2. 
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iiiiii • 


111 • 


1111 ■ 


1123 • 


13 • 


11 • 


ii 


00 • 


A A 
44 • 


OQ 
ZO 


ii • 


iz 


mm • 


111 ■ 


1111 • 


1111 • 


23 • 


13 • 


11 


• 11 • 


11 • 


11 


33 • 


44 


mm • 


111 • 


1111 • 


1111 • 


11 • 


23 • 


13 


• 11 • 


11 • 


11 


11 • 


11 


mm • 


111 • 


1111 ■ 


1111 • 


11 • 


11 • 


23 


• 13 • 


11 • 


11 


11 • 


11 


mm • 


111 • 


1111 • 


1111 • 


11 • 


11 • 


11 


■ 23 • 


13 • 


11 


11 • 


11 


mm • 


111 • 


1111 • 


1111 • 


11 • 


11 • 


11 


■ 11 • 


23 • 


13 


11 • 


11 


mm • 


111 • 


1111 ■ 


1111 • 


11 • 


11 • 


11 


■ 11 • 


11 • 


23 


13 • 


11 


mm • 


111 • 


1111 ■ 


1111 • 


11 • 


11 • 


11 


■ 11 • 


11 • 


11 


23 • 


13 


Here, • represents (8). 























Remark 3.3. Suppose pj = ui. for 1 < j < fc in a state p = pi 
state Too(p) = p'l^Si ■ ■ • ® p'^, p'j — ui. is valid for 1 < j < fc + 1. 



• • ®Pl- Then in the 
We postpone the inverse scattering formahsm of the dynamics to Theorem [ 



3.2. Counting particles and anti-particles. RecaU that a(C) is defined in (j2.24p . In 

our present context, it is the number of ah the particles and anti-particles within a box 
specified hy C, & Bi, where the term 2C,i means that a bound pair is regarded as a pair 
of a particle and an anti-particle (whose color is unspecified). The symbol a means all 
kinds of (anti-) particles. 

Let p = pi ® ■ ■ ■ ® PL be a state and write its time evolution as 

TUpi^- 



>Pl} 



P\ 



where p] G Bi^ . We write pj = P° = (Cj,i, • ■ • , Ci,n, Cj;„, • • • , C^.i) G Bi^ . For any elements 



ai, . . . , Or of {2, 3, . . . , n, n, 
(3.2) 



..,2,1}, we define the counting function 



Pai.....,aAp) 

t>l j = l 

where Cj-g = Cj,3 7 etc. The dependence on ai,...,ar enters the first term only. The 
indices in Pai,....a,- will always be arranged in the order 2, 3, . . . , n, n, . . . , 3, 2, 1. The 
second term is finite due to Remark 13.31 In fact the double sum may well be restricted 
to ^fj^i X]j=t+i where the nonzero contributions are contained. This region is depicted 



as the S W quadrant of the time evolution patterns like Example 13.11 and 13.2 



Pi P2 



PL 




(3.3) 

The first term in p.2p is the number of (anti-)particles with colors oi, . . . , contained 
in the top row which is the state p itself. The second term counts all kinds of particles 
and anti-particles in the hatched domain in p.3p H By the definition it follows that 

(3-4) Puip) = P2,...,«,n....,2a,T(7'oo(p))- 



^ More precisely, it should be hatched in a staircase shape. 
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Given a state p = pi ® ■ • • p/,, we write 
(3.5) P[k]=Pi(^---(^Pk il<k<L). 

Example 3.4. Let p be the state in the first fine in Example 13.21 Then, the counting 
function Pai,....ar{P[k]) 1 < fc < 9 takes the following values. (The middle column 
shows g such that pai,...,a^ = Pg in p.7|) - p.l0|) .) 



ai, . . . , Or 


9 


k = 1 


2 


3 


4 


5 


6 


7 


8 


9 
























23443211 


vo 


6 


11 


21 


32 


39 


46 


53 


60 


67 


2344321 


Vl 


5 


10 


19 


30 


37 


44 


51 


58 


65 


23443T 


V2 


4 


9 


17 


28 


35 


42 


49 


56 


63 


234T 




3 


8 


15 


24 


31 


38 


45 


52 


59 


2341 




2 


6 


13 


24 


31 


38 


45 


52 


59 


2T 


W2 — V2 


2 


5 


12 


20 


27 


34 


41 


48 


55 


T 


Wl — Vl 


1 


3 


9 


17 


24 


31 


38 


45 


52 





'^0 





2 


7 


15 


22 


29 


36 


43 


50 



We set Pai,...,ar(P[o]) = foi' ^-ny 0,1, ■ ■ ■ , cLt- Consider the difference P2i{P[k\) ~ Pr(P[fc]) 
for example. By the definition (j3.2p . it is the number of color 2 particles contained in 
Thus we have 

tt(color 2 particles) in pk {— Cfc,2) 
= P2-iP[k]) - PiiPlk]) - P2T{P[k~i\) + Pr(p[fc-i]). 

This is an example of the relations that reproduces a local variable from non-local count- 
ing functions. Given Z^, the set of counting functions that are necessary and sufficient to 
completely reproduce the local state pk G B/^ is not unique. However there is a choice 
that is linked with the generalized energies in Section 12.31 By using the function 7g in 
(I22Sl)-(E2ni), we set 



(3.6) Pa{p)^Y.^M + Y.Y. <P'. 



L L 

t\ 
j' 

j=l t>l j=l 



foi g = Va {0 < a < n — 1) , Wa ~ Va {I < a < n — 2), v!^_i and Vq^ . Although the 
notations pg here and Pai,...,ar in (|3.2p are somewhat confusing, we dare to use the both 
in the sequel supposing the resemblance is not too serious. Then p.6p is explicitly given 
as follows: 

(3.7) p,„(p) = P2,...,„,w,...— ,t(p) (0<a<n-2), 

(3.8) Pv^^^{p) = P2,3.....n-l.n,T(p). (p) = P2,3,...,n-l,rr,Tb)' 

(3.9) p^^_„^(p) = P2,3,...,aj{p) (l<a<'^-2), 

(3.10) Pvl^{p) = P${p)- 

The last one is subsidiary in that p^,"! (p) — p^i-viip) — Pvq(jp) + Pvi{p) holds reflecting 
(j2.12|) . One may also additionally introduce 

L L 

Pwa (p) = Pwa-va {p) + Pva (p) = X! fe) + 2 X! X! 

j=i t>i j=i 

for 1 < a < n - 2. See (|O0|) . For d[^\ the counting functions (|3J | - (f3?T0l) are precisely 
those listed in Example [ 
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Theorem 3.5. Forpjj.] in i3.5]) . setSpg = Pg{P[k]) ^ Pg{Plk-i])- The counting functions 
Jg. 9| )-( rO)) reproduce the local state pk = (Ci; • • • i Cn? Cni • ■ • 7 Ci) £ Bi^^ by 

Cl = h — Spvo + SPwi-vi , 

Ca = - <^P-ii.„_i-i,„_i (2<a<n-2), 

C„ = max(i5p„„_^ - (5p„._^,0), 
C„ = max((5p„;_^ - (5p„„_i,0), 

Ca = (^/ODa-i - '^Pt'a (1 < a < n- - 2). 
Proof. Straightforward by using ([S^), ([SJll-dSll) and niin(C„,C„) =0. □ 

4. Main result 
4.1. Counting functions and generalized energies. 

Theorem 4.1. For any state p G -B/j • • • Cs) i?;^ , the counting functions and the gener- 
alized energies i2.32\) coincide, namely, 

(4.1) £g{p)^Pgip) 

for g ^ Va {0 < a < n — 1) , Wa — Va {1 < a < n — 2) , and Uq^ . 

Here, should be understood as Sw^ip) ~ ^va{p)i the same convention is 

assumed for Iu,^-Va{p) in the sequel. Of course f^^(p) = PwSp) follows as a corollary. 
The (7's in Theorem 14.11 are the same as those considered in Lemma [2.41 (ii) and p. 61) . 
By substituting p.7p - p.l0p into (|4.ip . the theorem may be rephrased as 

£vAp) = P2,...,n,n,...,-^,iiP) (0 < a < n - 2), 

£v„^^ip) = P2,3....,n-l.n,Tb): (P) = ^2,3, . . .,n-l,n,T(p) ' 

Sw^-Va{p) = P2,3,....a,T(p) (l<a<"--2), 

For the proof we need one more Lemma. 

Lemma 4.2. Let p ^ pi® ■ ■ ■ ®pl € -B/i ® ■ • ■ (8) -B/j^ . For those g 's in Theorem \4.1\ the 

following equality is valid: 

L 

(4.2) Eg{p) - SgiT^ip)) =Y.g{i^'^ ^Pi), 

i=l 

where S i?oo is defined by Uao ®pi®---® Pi-i ~p'i®---® p-_j (g) ^^'^ . 

Proof. The following proof simplifies the one for Proposition 4.6 in Ref. [IT] in that the 
assumption ^i > • • • > Zl is not needed. We illustrate it for i = 3. Set p = pi ® P2 ® P^ 
and Too{p) = p'l ®-P3. Then, Proposition 12.61 tells that £g{Tao{p)) is the sum of g at 
all • in the following diagram. 

Uoo Uoo Pi P2 P3 
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On the other hand, the right hand side of (|4.2p is equal to the sum of g at all o. Thus 
from Lemma 12.41 (ii), we find 

L 

^g{Too{p))+^9{i^'^ ®Pi) = ^ Ig{Uoo®Uoo®Pl®---®Pj) 

i=l l<j<3 

= 2^ /g(Uoo «)pi ® • • • = £g[p), 

1<J<3 

completing the proof. □ 

Proof of Theorem \4.1\ From Remark l3.3[ T^{p) = w;^ ® • • -^ui^ for t > L. For such a 
state = and Pg = hold due to (|2.23p and p.7p - p.l0p . respectively. Thus it suffices 
to show 

Sg{p) - £g{Too{p)) = Pg{p) " Pg{Too{p))- 

By applying (j4.2p and (13.61) to the left and the right hand sides, respectively, this becomes 

L L 

®P^) = H^lsiP^) + a(pD - lg{p[)), 
1=1 1=1 

where we have set Too{p) — p'l <E) ■ ■ ■ ^ p'j^- (This was denoted hy p\ (g) ■ ■ ■ ® p\ in p. 61) .) 
From the definition of ^('^ in (|4?2|) . we have ^^^^ ~ p - (K) C*-'^^'- Therefore Lemma [2^ 
tells that = lg{Pi) + ci(Pi) ~ 7g(-Pi) holds for each i, finishing the proof. □ 

4.2. ^-transformed correspondence. Let us give an analogous result on 5 = (1 < 
a < n — 2) which is not included in Theorem 14. II Our presentation in this subsection is 
brief since the essential features are the same as the previous case. To state the result, 
let us introduce a *-transformed generalized energy and a ^-transformed Dn^ cellular 
automaton. (See l|2.9p and (|2.10l) for the original definition of *.) 

Let u* = (0, . . . , 0, /) e i?/. See (|2.17l) . The ^-transformed generahzed energy £*, {pl® 
■ ■ ■®Pi) of an element pl® ■ ■ ■ ®Pi G Bi^ ■ • • Bi^ is the sum of w* for all the vertices 
in the following diagram (L = 3 example): 

P3 P2 Pi u*^ 

£*.{P:i®P2®Pi) = 




Compare this with Proposition 12.61 One can show that E*, is well defined and R- 
invariant. 

The ^-transformed d\1^ cellular automaton is the dynamical system on Bi^ ®---® Bi^ 
endowed with the commuting time evolutions Ti{l > 1) defined by (g) • • • i^) Pi CE> u* — 
^ (g) Ti*{pl ® • ■ • g) Pi)- {£, G Bi is determined by this relation.) is well defined. 
Moreover, under the time evolution p'^ g) • • • ® p'l = T^{pL €5 • • • ig" Pi), the equality 
p'j = u*. is valid for 1 < j < fc-|- 1 if pj = u*i. for f < j < A:, which is parallel with Remark 
For ( E Bi, introduce the charge conjugation of (|2.24l) - (|2.25p by 

a*(C) a(C*)-2Ci+C2 + --- + Cn + C„ + --- + C2, 

7:.(C) 7..(C) = Ci+Ca+i + --- + C« + C„ + --- + C2 (l<a<n-2). 
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Writing the time evolutions of a state p — (E) ■ ■ ■ Pi € Bi^^ ® ■ ■ ■ ® Bi^ as {T^Y{p) — 
Pi^ <E) ■ ■ ■ <E p\, we define the counting function (f < a < n — 2): 

PL ■■■ P2 Pi 

L L 

j=i t>i j=i 



The counting is done by 7*. for the top row and by a* for the SE quadrant generated by 
beneath it. Thanks to the commutativity of * and the combinatorial R (Prop. 4. 4 in 
Ref. [IS]), Theorem 14.11 implies the following: 

Proposition 4.3. For any state p e Bi^^ • • • ig) Bi^ , the following equality is valid: 

£:,{p) ^ P:.{p) (l<a<n-2). 

This completes our interpretation of the (^-transformed) generalized energies associ- 
ated with all the generalized local energies in terms of the (^-transformed) Dn'^ cellular 
automaton. 

5. Connection with combinatorial Bethe ansatz 

Combinatorial Bethe ansatz was initiated by Kerov-Kirillov-Reshetikhin (KKR) [121 
[13] to establish a fermionic formula of the Kostka-Foulkes polynomials with the invention 
of rigged configurations and the KKR bijection. Their fermionic formula generalized 
Bethe's formula [4 for some simplest Kostka numbers that originates in the completeness 
issue. 

Rigged configurations are combinatorial analogue of solutions to the Bethe equations. 
The KKR bijection maps them to the combinatorial analogue of Bethe vectors which may 
be viewed as elements of (a subset of) Bi-^ ® • ■ ■ ® Bi^ . The combinatorial Bethe ansatz 
has flourished in the fermionic formulas for general affine Lie algebras [Sj \19\ [20l [21] , the 
solution of the initial value problem of integrable An^ cellular automata by the inverse 
scattering method and a connection with the classical soliton theory |5I via ultra- 
discrete tau functions |17| and so forth. Our aim in this section is to present an inverse 
scattering formalism of the Z?!^"* cellular automaton and to conjecture explicit formulas 
for some generalized energies in the form of ultradiscrete tau functions associated with 
the Dn'' rigged configurations. 

5.1. Inverse scattering formalism. Set V+ — {p € Bi^ ® ■ ■ ■ ® Bi^ \ iip = for i — 
1,2,..., n}. A state belonging to V+ is called highest. It is known that there is a bijection 
between Vj^ and the set of rigged configurations [201 121] . Consider a set 

(5.1) 5 = {(a„j„r,) e {l,2,...,n} x Z>i x Z>o \ i = 1,2, . . . ,N}, 

where iV > is arbitrary and each triplet s = (a, j, r) called string possesses color, length 
and rigging which will be denoted by cl(s) = a,lg(s) = j and rg(s) = r, respectively]! 
5 is a rigged configuration if rg(s) < p[g^(^g^'' is satisfied for all s E S. Here pj"^ = 
^a,i Ylk=i niin(j, lk)-J2tes ^a,ci(t) mhi(j, \g{t)), where {Ca,b)i<a,b<n is the Cartan matrix 
of Dn- Note that PigJ^j')^^ ^ has to be satisfied for all s d S, which imposes a stringent 
condition on the set {(a^, ji) \ i = 1, . . . , N}. Set RC — {S : rigged configuration}. 

■^Colors l,2,...,n of strings in rigged configurations should not bo confused with colors 
1,2, ... ,n,n, ... ,2,1 of (anti)-particles. 
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Theorem 5.1 (Refs. ^ and [H]). There is a bijection $ : P+ RC. 

An explicit algorithm to determine the image of is known. It is a Dn ^ analogue 
of the Kerov-Kirillov-Reshetikhin bijection [T3] for An \ which plays a central role in the 
combinatorial Bethe ansatz. Our convention here is the one adopted in Ref. ^1^. 

For any highest state p S V+, its time evolution Ti{p) is again highest. Thus T/ induces 
a time evolution on RC via $. Let L be large enough and assume the situation in (j3.1l) . 
Then we have 

Theorem 5.2. ^{Tiip)) = fi<^{p) holds, wherefi : {{a^,j^,ri)} ^ {{a^, j^,r^+5aiSm.l■a{j^, 
is a linear flow on rigged configurations. 

Proof. The proof uses Theorem 8.6 of Ref. [21] and is similar to that of Proposition 2.6 
in Ref. [H] for type aI^\ □ 

Thus the composition "I>^^ oT; o $ linearizes the original time evolution T; and solves 
the initial value problem in the Dn^ cellular automaton by the inverse scattering method. 
See Ref. [14] for an analogous result for An \ 

5.2. Conjecture on ultradiscrete tau functions. For a rigged configuration S (|5.ip . 
let T C be a (possibly empty) subset of S. In general, T is no longer a rigged 
configuration. We introduce the piecewise linear functions {0 < k < L and < d < n) 

= ^ 51 ^-i(«).-i(t)™n(lg(s),lg(t)) + ^rg(s), 

s.teT sST 
k 

cf(T) = c{T)-J2 E niin(?„lg(.))+ ^ lg(.) 

1=1 s(ET,cl(s) = l seT,cl{s)=d 

By the definition, the last term in c'^\t) is when d = 0, and the relation 

(5-2) C^k^ (T) ^ 4"^ (r) |rg(s)^rg(s)+lg(s)A-,i(,,,a 

holds. Obviously we have c(0) — c"^\%) = 0. On the other hand, c{S) is known as the 
(co)charge of the rigged configuration 5 [OlISl l^l^ . 

We define a Z>o-valued piecewise linear function on S as follows: 

(5.3) T('^)(5)--min(4'^)(T)) (0 < A: < L, < d < n). 

For S in (jS.ip . the minimum extends over 2^ candidates and reminds us of the structure 
of tau functions in the theory of solitons [8 . In fact, for type An \ analogous functions 
have been identified [TTl as ultradiscretization of the tau functions in KP hierarchy. 
Although such an origin is yet to be clarified, we call (|5.3p ultradiscrete tau function. 
Guided by the results in An'^ and supported by computer experiments, we propose 

Conjecture 5.3. For any highest state p G V^, let S — ^{p) he the corresponding rigged 
configuration. Then, the following equalities hold for p\j^^ \3. 5]) with < k < L. 



(5.4) 




= ^vo{P[k]), 


(5.5) 




= Sv^^iP[k]), 


(5.6) 




= Svl_^{P[k]), 


(5.7) 




= ^Vr,-l{Plk]), 


(5.8) 




= {P[k] ) - {P[k] ) + ifiO {P[k] ) 
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In (jS.Sp . ^Q(j>[k]) is the standard notation in crystal theory meaning niax{j > 
/oP[fc] ^ 0}. By using (EH), dSH]) with a = 0, (jXTUI) . and Theorem O one can 

show that (|5.4p and (|5.5p are equivalent. 

The algorithm [20, 21 for seems valid not only for highest but arbitrary states 
if one allows negative rigging. With such a generalization, we expect that Theorem 15.21 
and Conjecture 15.31 hold for any state, which was indeed the case for type A^^^ [H]. 

Example 5.4. For the initial state p in Example 13.21 its rigged configuration is S' = 

= {(1,8, -2), (1,6,0), (1,2, -1), (1, 1, -1), (2,8,0), (2,6, -1), (2,2, -1), (3,8, -3), (4,8, 
The ultradiscrete tan function t^^\s) and (pa{p[k]) take the following values. 





k = 1 


2 


3 


4 


5 


6 


7 


8 


9 




6 


11 


21 


32 


39 


46 


53 


60 


67 


ri'\s) 





2 


7 


15 


22 


29 


36 


43 


50 




1 


3 


9 


16 


23 


30 


37 


44 


51 


ri'\s) 


2 


6 


13 


24 


31 


38 


45 


52 


59 


rtHs) 


3 


8 


15 


24 


31 


38 


45 


52 


59 


fo{p\k]) 


1 





1 





















Comparing this with Example 13.41 one can check Conjecture 15.31 

Still many generalized energies in previous sections await formulas as in Conjecture 
15.31 to be discovered. They are ultradiscrete analogue of the so called X = M conjecture 
[3 [19] in the sense that the generalized energies from crystal theory acquire explicit 
formulas of a fermionic nature originating in the combinatorial Bethe ansatz. Such 
results combined with Theorems 13.51 and 14.11 will lead to a piecewise linear formula for 
the bijection $^ as was done for A^^^ [T^. 

Let us end by raising a closely related question as another future problem. Let V+ (A) 
denote the subset of P+ having the prescribed weight A. Then, does the generating 
function 

X,(A)= E l'''-'^ 

pGP+(A) 

of the generalized energy admit a fermionic formula like Refs. ITS' and [S]? 
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